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Abstract. We show that the direct Urea process of neutrino emission is allowed in two possible phases of non- 
spherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. 
The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear 
structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the 
neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2-3 orders of magnitude 
higher than the emissivity of the modified Urea process in the stellar core. The process may noticeably accelerate 
the cooling of low-mass neutron stars. 
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1. Introduction 



It is well known (Lattimer et al. I1991|) that direct Urea 
process produces the most powerful neutrino emission in 
the inner cores of neutron stars (NSs) . The simplest direct 
Urea process in a dense degenerate matter, composed of 
neutrons (n) with an admixture of protons and electrons 
(p and e), consists of two successive reactions (direct and 
inverse ones): 



P 



P 



n+^e, 



(1) 



where v e and v c are the electron neutrino and antineu- 
trino, respectively. The process is allowed by momentum 
conservation if p Fn < p Fp + p Fc , and forbidden other- 
wise. Here, p Fn ,p Fp , and p Fc are the Fermi-momenta of 
n, p, and e. The neutrino momentum p v ~ k^T is much 
smaller than these Fermi momenta, and it can be ne- 
glected in momentum conservation (T is the internal NS 
temperature, and fee is the Boltzmann constant). It turns 
out that the direct Urea process is allowed only at suffi- 
ciently high densities (typically, a few times of po, where 
po = 2.8 x 10 14 g cm -3 is the nuclear matter density at 
saturation) for those model equations of state of dense 
matter which have a large symmetry energy (rather high 
fractions of protons and electrons) to satisfy the inequality 

PFn<PFp+PFc- 

Thus, the direct Urea process is forbidden in the cores 
of low-mass NSs. The main neutrino emission from these 
(nonsuperfluid) cores is produced by the modified Urea 
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process (nN— >peND e , peN—>nNi/ e , where N=n or p is a 
nucleon-spectator required for momentum conservation). 
The modified Urea process is 6-7 orders of magnitude 
weaker than the direct Urea process. This greatly reduces 
the neutrino emission of low-mass NSs in comparison with 
the emission of massive NSs, where the direct Urea is open. 

In this paper we analyze the possibility to open the 
direct Urea process in the inner NS crust. The main idea 
is that a momentum excess may be absorbed by a lattice 
of nonuniform nuclear structures in the crust. Specifically, 
we consider a model of nonspherical nuclear structures 
which appear in the density range from p sa 10 14 g cm~ 3 
to the crust-core interface (p cc sw po/2) if one employs 
some models of nucleon-nucleon interaction (Ravenhall et 
al.GHESI Pethick & R.avenhall HM)|l . The theory predicts a 
sequence of phase transitions with increasing p within this 
density range: from familiar spherical nuclei to cylindrical 
nuclear structures, from cylinders to slab-like structures, 
from slabs to inverted cylinders, then to inverted spheres 
and, finally (at p = p cc ), to a uniform nuclear matter 
in the core. The shell of nonspherical nuclei in the crust, 
sometimes called the NS mantle, is thin (not thicker than 
several hundred meters) but contains a noticeable frac- 
tion of the crust mass. We will consider two last phases 
- the inverted cylinders and the inverted spheres, where 
free protons appear (in addition to free neutrons in the 
inner crust), but periodic nuclear structures are still not 
dissolved into the uniform nuclear matter. The periodic 
structures modulate motion of neutrons and protons (in- 
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ducing Bloch states) and open direct Urea process in a 
neutron-star mantle. 

2. Periodic potential and nucleon wave functions 

In order to calculate the neutrino emissivity of direct Urea 
process in the NS mantle we need wave functions of neu- 
trons and protons in a periodic nuclear potential. For un- 
derstanding the main features of the problem, we adopt 
a simplified Thomas-Fermi approximation. The nuclear 
structures will be described using the results of Oyamatsu 
( 1993). Specifically, we will employ his model I for the en- 
ergy density functional. The potential energy of neutrons 
or protons (j = n or p) can be calculated as 



where we have used the identity 



cell 



(6) 



Vj(r) 



,5/3 



— — n 5/3 
2m n ? 



(2) 



where eo(n n ,n p ) is the energy density of the uniform nu- 
clear matter, rrij is the nucleon mass, and rij(r) is a local 
number density of nucleon species j, which depends on dis- 
tance r from the nucleus center inside a Wigner-Seitz cell 
[Eq. (4.8) of Oyamatsu[T993 . Equation (0) neglects small 
gradient corrections. For simplicity, we restrict ourselves 
to scalar (spin-independent) nuclear potentials. 

Using the Schrodinger perturbation theory we can ex- 
press the Bloch wave function of a nucleon (n or p) in the 
periodical potential Vj(r) as: 



Now in Eq. (|SJ) we can integrate over the simplified 
Wigner-Seitz cell (a cylinder or a sphere for the phases 
of inverted cylinders and inverted spheres, respectively). 

Unless the contrary is indicated, we will use the units, 
where ft = c = AsB=V = l. Notice, that the pertur- 
bation expansion J3J) fails (becomes singular) at Bragg's 
diffraction points (at \p\ = \p + q\) which signals the spe- 
cial importance of the band structure near these points. 
However, we will see that nucleons with such "resonant" 
wave functions do not contribute to the neutrino process 
of study. 

3. General formalism 

The emissivity Q of the direct Urea process in the NS 
mantle is calculated in the same way as in the NS core 
(Lattimer et al. 1991). Using the notations from the review 
article by Yakovlev et al. (|200T). we obtain: 



Q = 2 



(2tt) 3 



dW^f e„/ n (l-/p)(l-/ e ), 



(7) 



JQ 



En — E„> 



where q is an inverse lattice vector, p' = p + q, V is the 
normalization volume, Co = 1, and C q = Vj q /{E p — E p >) 
for q ^ 0. Furthermore, \s is a unit spinor (XsXs' = Sss'), 
s = ±1 is the sign of the nucleon spin projection onto 
the quantization axis, E p = p 2 /2m* is the unperturbed 
energy, pj is the momentum, to* is the effective mass at 
the Fermi surface, and Vj q is a Fourier-transform of the 
potential Vj(r): 

' ' V 3 {r)e- lqr dr, (4) 



where fj is the Fermi-Dirac distribution of particle species 
j(=n, p, e), e v is the neutrino energy, dW,-,/ is the dif- 
ferential probability of neutron decay [calculated with the 
wave functions 10)]. The overall factor 2 doubles the emis- 
sivity of the neutron decay to account for the contribution 
of the inverse reaction (assuming /3-equilibrium). After 
— pla ^y^C q e tp r , (3) standard simplifications (e.g., Yakovlev et al. I2(JU1|) and 
* q integration over orientations of the neutrino momentum, 

we get: 



dW, 



dp n 



(27T) 4 



(2tt) 3 (2tt) 12 

x 5(p' n -p' p - Pc ) 
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x \M }1 \ 2 At: e 2 v de v \\p Fi 

3=1 



m* dej dfl/, (8) 



cell 



Vjq = 

VceH 

Vceii being the elementary-cell volume. An exact calcula- 
tion of Vj q is difficult because the configuration of the ele- 
mentary cell is complicated. However, the calculation can 
be simplified due to the short-ranged character of nuclear 
forces. Consequently, the nucleon potential is almost con- 
stant near the cell boundaries, where it can be replaced by 
the boundary value Voo ■ Introducing a simplified Wigner- 
Seitz cell of equivalent volume as defined by Oyamatsu 
fiMfy . p. 434, we find: 



where ej is the energy of particles species j, dflj is the 
solid angle element in the direction of Pj, \pj\ = pp j (i.e., 
nonperturbed particle momenta pj in the momentum- 
conserving delta function are placed at Fermi surfaces), 
and m* = jU e {p-o being the electron chemical potential). 
Finally, 



\M }1 \ 2 = 2 G 2 cos 2 0c \C q J 2 \C qp \ 2 (j 2 + 3g 2 A ) 



(9) 
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(^■(O-VooJe-^dr, 



(5) 



is the squared matrix element, summed over particle spins 
and averaged over directions of neutrino momentum. Here, 
Gf = 1.436 x 1CU 49 erg cm 3 is the Fermi weak interaction 
constant, #c is the Cabibbo angle, sin#c = 0.231, fy = 1 
is the vector interaction constant, and <?a = 1.26 is the 
Gamow- Teller axial- vector interaction constant. 

The leading term in the sum © over inverse lat- 
tice vectors corresponds to q n = q p — 0. However, in 
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this term p' n = p n , p p = p p (see Eq. (0 and the dis- 
cussion afterwards), and the neutrino emission is forbid- 
den (p^ > p Fp + p Fc ) by momentum conservation. Thus, 
the main contribution into the emissivity Q comes from 
smaller terms with either q n = or q p = 0. The terms 
with q n 7^ and q p ^ are even smaller and can be ne- 
glected. The retained terms are constructed in such a way 
that the momentum-conserving delta function excludes 
the "dangerous" Bragg diffraction points. For instance, 
one can easily show that for q n = the Bragg diffraction 
condition \p p \ — |p p + g p | is incompatible with momentum 
conservation in Eq. JHJ. 

Inserting Eq. © into Eq. Q and integrating over par- 
ticle energies and propagation directions, we get: 



Q(T,p) = Q Q (T,p)R(p), (10) 
Qo(T, p) = igl G% cos 2 9 C (/* + 3 5 2 ) m> p TO e * T 6 



, i n , \ ^ 3 m*m* erg , s 

4.0xl0 27 ( — I -^T 9 6 -f-, (11) 



R(p) = E E 
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where n c is the electron number density, no = 0.16 fm 3 
is the nucleon number density in saturated nuclear mat- 
ter, T 9 = T/10 9 K, aj = q/ PF] , Df n = max(-l, D 3 -_), 
Df** = min(l, D j+ ), 



F(x) 
D n ± 
D P ± 



s/TTx + l 



VT+x- l 



-m 

2 



(Pfp^Pfc) 2 -^'^ 
2 Pf„9 

(Ppn±PFc) 2 -PFp 



2 PFp9 



with p Fn - p Fp - p Fo < q < p Fn + p Fp + p Fc [and F(x) = 
otherwise] . In Eq. l|10fl Qo is the direct Urea emissivity in 
the uniform nuclear matter (disregarding the momentum- 
conservation constraint) and R can be called its reduction 
factor which describes the weakening of the process in the 
NS mantle. This factor depends on p, but not on T. The 
emissivity Q is oc T 6 , just as for the well-known direct 
Urea process in the NS core. 

For calculating R from Eq. Ijl2(l we need the parti- 
cle Fermi-momenta. Because our analysis is simplified, we 
have used two models. First, we have defined the Fermi- 
momenta by: 



PFn 



(37r 2 n„) 1/3 , p Fp = p Fc = (3Tr 2 n p ) 



1/3 



(13) 



where n n and n p are the nucleon number densities aver- 
aged over the Wigner-Seitz cell [see Eq. (4.8) and Table 
6 of O vamatsu 11993) . and p Fp = p Fc due to electric neu- 
trality. Second, we have adopted electric neutrality and 



/3-equilibrium, and determined the Fermi-momenta from 
the equations: 



P 



Fn 



Pfc-Pfp- 2mn 



Ko = 
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Fp 



2m r 



V p0 +p Fe , 



(14) 



where Vjo is the central (q — 0) Fourier harmonics of the 
nucleon potential Vj (r) . 

The values of R(p) calculated from Eq. I|12|) appear to 
be qualitatively the same for both models. In our calcula- 
tions we have set m* = mj , but variations of m* within 
reasonable limits do not qualitatively change R(/o). The 
results based on the model (|14fl can be fitted by 



R(p) w R 2 + (Ri-R 2 )(l-x) 2 x<l, 
R(p) w R 2 /x 7 x>l, 



(15) 



where Ri = 6x 10 -5 , R 2 = 10~ 5 , x = (rib — ni) / (nz — ni) , 
n b = p/m n is the baryon number density, while n\ and 
n 2 are the baryon number densities at the outer and in- 
ner boundaries of the inverted-cylinder phase. Thus, the 
density parameter x varies from to 1 in the layer of 
inverted cylinders, and we have x > 1 in the layer of in- 
verted spheres. In the Oyamatsu model, which we em- 
ploy, m = 0.08274 fm" 3 , n 2 = 0.08537 fm~ 3 , and the 
layer of inverted spheres extends to na = 0.08605 fm~ 3 
(to x = 4.868). 

The number N of inverted lattice vectors, which con- 
tribute into R, is large: N ~ 200 for the phase of inverted 
cylinders, and N is up to ~ 2800 for the phase of in- 
verted spheres. The proton contribution (q p 0) into R 
is approximately three times larger than the neutron one 
(Qn + 0). 

As seen from Eq. (|15|l . R ~ 10 5 . Thus, the emissivity 
of the direct Urea process in the mantle is about 5 orders of 
magnitude weaker than in the inner NS core. Nevertheless, 
as will be shown in the next section, the direct Urea in the 
mantle can affect the cooling of NSs. 

4. The effect on the cooling of low-mass NSs 

We will focus on sufficiently low-mass NSs, with the for- 
bidden direct Urea process in the inner cores (Sect. [JJ. 
The cooling behaviour of these stars is described, for in- 
stance, by Potekhin et al. (2003). The neutrino luminosity 
L v of low-mass stars is not too high. Thus, the additional 
neutrino emission from the mantle may be pronounced at 
the neutrino cooling stage (when the stellar age t < 10 5 
yrs). 

An order-of-magnitude estimate gives 



L v ~ 4 7T i& h Q + (4vr/3) i? 3 c Q c 



(16) 



where R cc is the core radius, h is the mantle width, 
and Q cor c is the mean core neutrino emissivity. Adopting 
Rcc = 10 km and h = 100 m, we obtain that the direct 
Urea in the mantle may affect the cooling if Q J> 30 Q CO re, 
which is quite possible. For instance, this inequality may 
hold for non-superfluid low-mass NSs. The main neutrino 
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Fig. 1. Critical temperature T c versus p for various mod- 
els of neutron and proton pairing in the NS core and crust. 
Models C93 and S: singlet-state neutron pairing in the 
crust; models lp and 2p: singlet-state proton pairing in 
the core; model lnt: weak triplet-state neutron pairing in 
the core. The crust-core interface and the threshold den- 
sity of the direct Urea process in the core are shown by 
vertical dot-and-dashed lines. 



emission in the cores of these NSs is produced by the modi- 
fied Urea process, which is 6-7 orders of magnitude weaker 
than the direct Urea process in the cores of massive NSs, 
while the direct Urea process in the mantle is five orders 
of magnitude weaker than in the cores of massive NSs. 

Let us illustrate these statements by more elaborated 
calculations. For this purpose we employ the equation of 
state of Prakash et al. (1988 ) in the stellar cores (the ver- 
sion with the compression modulus of the saturated sym- 
metric nuclear matter K — 240 MeV and with model I 
for the symmetry energy). The NS models based on this 
equation of state are described, e.g., by Yakovlev et al. 
(2001}. The equation of state opens direct Urea process in 
the NS core at p > 7.85 x 10 14 g cm~ 3 (which is possible 
in NSs with M > 1.358 Mq). The most massive stable 
NS has the central density p c = 2.578 x 10 15 g cm -3 and 
the mass M = 1.977 Mq. A typical low-mass NS with 
M = 1.35 M has p c = 7.79 x 10 14 g cm" 3 and R = 13.0 
km. 

We have simulated the NS cooling with our fully rel- 
ativistic cooling code (Gnedin et al. 1200 l|l . It calculates 
the cooling curves: the NS surface temperature T s °°, as 
detected by a distant observer, versus the NS age t. We 
have updated the code by incorporating the neutrino emis- 
sion from the direct Urea process in the NS mantle. We 
have considered nonsuperfluid NSs and NSs with neutron 
and proton superfluidity of the internal layers. We have 
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Fig. 2. Total neutrino emissivity versus density at T = 
3 x 10 8 K. The crust-core interface and the direct Urea 
threshold in the core are marked by dot-and-dashed lines. 
The peak near the core-crust interface shows the direct 
Urea process in the mantle (presented also in the insert). 
Dotted line: no nucleon superfluidity (nosfl); solid line: 
model lp of proton superfluidity; dashed line: model lp 
of proton superfluidity in the core and model S of neu- 
tron superfluidity in the crust (see the text). The neutron 
pairing lnt is too weak to appear and affect the neutrino 
emission at the given T . 



taken into account (i) a possible singlet-state pairing of 
free neutrons in the crust and the outermost part of the 
core, (it) a triplet-state pairing of neutrons in the core, 
and (Hi) a singlet-state pairing of protons in the crust. 

Microscopic theories of nucleon superfluidity of dense 
matter give very model dependent density profiles of su- 
perfluid critical temperatures of nucleons, T c (p) (e.g., 
Lombardo & Schulze I2001f) . Thus, we have considered 
several superfluid models (Fig. 1) available in the liter- 
ature: one phenomenological model lp of strong singlet- 
state proton pairing, and one phenomenological model 
lnt of weak triplet-state neutron pairing in the NS core 
(Kaminker et al. l2002[> : models S (Schulze et al.[T996) and 
C93 (Chen et al. 11993(1 of singlet-state neutron pairing 
in the crust. We have also proposed the additional phe- 
nomenological model 2p of proton pairing in the core (not 
to be confused with model 2p in Kaminker et al. [2002!). 
This pairing is rather weak at the core-crust interface but 
becomes much stronger at higher p. 

The effects of superfluidity on the neutrino emission 
and heat capacity of the matter have been included in 
the standard way (Yakovlev et al. [1999 ) . It is well known 
that the nucleon superfluidity reduces the traditional neu- 
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Fig. 3. Theoretical cooling curves compared with obser- 
vations of two NSs. All the lines but the short-dashed one 
are for 1.35 Mq NS (direct Urea in the core forbidden). 
Any pair of lines of the same type refers to one super- 
fluidity model with the direct Urea process in the mantle 
on (lower curve) or off (upper curve). Dotted curves: non- 
superfluid NSs (nosfl); solid curves: lp and lnt superflu- 
idities in the core; long dashes: lp and lnt superfluidity 
in the core and C93 superfluidity of neutrons in the crust; 
dot-and-dash lines: 2p and lnt superfluidities in the core. 
The short dashed line is for a nonsuperfluid 1.5 Mq NS 
with the open direct Urea process in the core. 



trino mechanisms but opens the neutrino emission due to 
Cooper pairing of nucleons. We assume that free protons 
and free neutrons in the mantle have the same T c (p) as 
protons in the core and free neutrons in the ordinary crust 
of spherical nuclei. We have adopted the same form of the 
superfluid reduction factor of the direct Urea process in 
the mantle as in the inner core. 

Figure 2 shows the density profile of the total neutrino 
emissivity at T — 3 x 10 s K for nonsuperfluid and super- 
fluid NSs, while Fig. 3 presents some cooling curves. All 
the curves, but one short-dashed curve, are calculated for a 
1.35 Mq NS with the forbidden direct Urea process in the 
core. This is a typical example of a low-mass NS. Any pair 
of lines of the same type corresponds to the same super- 
fluidity model. Every upper line of the pair is calculated 
neglecting the direct Urea process in the mantle while its 
lower counterpart is calculated including this process. The 
short dashed curve is an example of the cooling of a more 
massive, 1.5 Mq, nonsuperfluid NS. Its central density is 
p c = 9.0 x 10 14 g cm -3 ; the direct Urea process is open in 
the inner stellar core, producing a strong neutrino emis- 
sion and a rapid cooling of the star. In this case the direct 
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Urea process in the mantle is entirely negligible: it cannot 
compete with its senior partner in the core. 

Along with the theoretical curves, in Fig. 3 we present 
the observational limits on the effective surface temper- 
atures and ages of two NSs, RX J0822-4300 and PSR 
B1055-52. Among all middle-aged isolated NSs, whose 
thermal emission has been observed and the effective sur- 
face temperature measured (the data are given, e.g., in 
Yakovlev et al. 2003), these are the NSs hottest for their 
ages. They can be interpreted as low-mass NSs (see, e.g., 
Yakovlev et al. RjTJOcSl and references therein) . 

RX J0822-4300 is a radio silent NS, a compact central 
object in Puppis A. Its effective temperature is taken from 
Zavlin et al. <jTM9"|l . Recently Pavlov lj2TjQ"3|> has kindly 
provided us with the updated value of the effective tem- 
perature of PSR B1055-52, T s °° « 7 x 10 5 K. However, he 
has not indicated errorbars. We introduce (Fig. 3), some- 
what arbitrarily, 10% uncertainties of T s °°. The ages of RX 
J0822-4300 and PSR B1055-52 are taken as described by 
Yakovlev et al. <|20"0"5|) . 

The dotted line in Fig. 2 shows the neutrino emissivity 
in a nonsuperfluid NS. The peak before the crust-core in- 
terface is produced by the direct Urea process in the man- 
tle (additionally presented in the insert). The emissivity 
at the peak maximum is about three orders of magnitude 
larger than in the NS core. The appropriate cooling curves 
are shown by the dotted lines in Fig. 3. The direct Urea 
process in the mantle noticeably accelerates the cooling at 
t ~ 10 5 yrs. 

Assuming the strong lp proton pairing and the weak 
lnt triplet-state neutron pairing in the core (and the man- 
tle) but no singlet-state neutron pairing in the crust, we 
obtain the emissivity profile plotted by the solid line in 
Fig. 2. The proton superfluidity fully suppresses the mod- 
ified Urea process in the outer core, before the direct Urea 
threshold, and partly suppresses the direct Urea process in 
the inner core (of massive NSs) beyond the threshold. The 
neutron pairing in the core is so weak that it does not ap- 
pear at T = 3x 10 s K in Fig. 2. The neutrino emission from 
the core of a low-mass NS becomes very slow, being mainly 
determined by the neutron-neutron bremsstrahlung pro- 
cess. This increases the surface temperatures of middle- 
aged NSs (the solid curves in Fig. 3), and enables one 
to interpret the observations of RX J0822-4300 and PSR 
B1055-52 (e.g., Kaminker et al. l2"07)2l Yakovlev et al. 
2003)). In this case the proton superfluidity partly sup- 
presses the direct Urea process in the mantle and reduces 
the difference of the cooling curves calculated with and 
without this process. Notice that the combination of the 
superfluidities lp and lnt is sufficient to explain the data 
on all isolated middle-aged NSs [whose surface temper- 
atures have been measured (estimated) from the obser- 
vations of their thermal radiation] by theoretical cooling 
curves of NSs with different masses (e.g., Kaminker et al. 
2002 ) . We sec that the inclusion of the direct Urea process 
in the mantle does not violate this interpretation. 

If we additionally switch on the singlet-state super- 
fluidity C93 in the crust, we will obtain the long-dashed 
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cooling curves in Fig. 3, which go slightly lower than the 
solid curves. This additional acceleration of the cooling is 
produced by the neutrino emission due to Cooper-pairing 
of neutrons in the crust. The effect is weak because the 
superfluidity C93 dies out long before the crust-core in- 
terface and produces the weak Cooper-pairing emission 
(Potekhin et al. 12003)1 . The relative importance of the di- 
rect Urea process in the mantle is even lower than in the 
absence of the superfluidity C93. 

If we add the singlet-state neutron superfluidity S in- 
stead of C93, the situation would be different in two re- 
spects. First, the superfluidity S completely switches off 
the direct Urea process in the mantle (see the dashed curve 
in Fig. 2; its peak in the core is produced by the neutrino 
emission due to the pairing S). Second, the superfluidity S 
occupies much larger fraction of the NS volume than the 
superfluidity C93, intensifying the neutrino emission due 
to the singlet-state pairing of neutrons. Now this emission 
noticeably lowers the cooling curve, complicating the in- 
terpretation of RX J0822-4300 and PSR B1055-52. We 
do not present the cooling curve for that unlikely case in 
Fig. 3; it is the same as given by Potekhin et al. (2003). 

Finally, let us employ the proton pairing 2p and the 
neutron pairing Int in the core (and neglect the singlet- 
state neutron pairing in the crust). The proton superfluid- 
ity 2p suppresses the neutrino emission in the main frac- 
tion of the core but not in the mantle. The effect of the 
direct Urea process in the mantle on the cooling of NSs 
of the age t ~ (1 4 — I0 5 ) yrs becomes most pronounced 
(the dot-and-dash curves in Fig. 3). With this process on 
(the lower dot-and-dash curve) the low-mass star will cool 
too fast, strongly complicating the interpretation of the 
observations of RX J0822-4300 and PSR Bf 055-52. 

Let us emphasize that the existence of the NS mantle 
(the layer of nonspherical atomic nuclei) is still a hypothe- 
sis. The theory predicts this layer only within some models 
of nucleon-nucleon interaction. Since the lower dot-and- 
dash curve in Fig. 3 strongly contradicts the observations, 
the underlying physical scenario becomes doubtful. This 
would imply, for instance, that the 2p proton superfluid- 
ity model is inadequate, or the neutrons in the mantle are 
strongly superfluid (switching off the direct Urea process), 
or the mantle does not exist at all. 

Our analysis of cooling low-mass NSs is illustrative and 
incomplete. The cooling of these NSs is actually affected 
by (?) NS superfluidity, (it) NS surface magnetic fields, 
and (Hi) possible surface layer of light (accreted) elements 
(as discussed in detail by Potekhin et al. 12003)1 . Our cal- 
culations indicate the existence of the fourth regulator, 
(iv) the presence of the NS mantle and the associated di- 
rect Urea process. As clear from the results of Potekhin et 
al. 1)2003)1 and our present results, all four regulators are 
of comparable strength and should be analyzed together. 
This many-parametric analysis is beyond the scope of the 
present article. 



5. Conclusions 

We have calculated the neutrino emissivity of the new neu- 
trino mechanism - the direct Urea process in a neutron- 
star mantle, a thin layer of nonspherical nuclei (Ravenhall 
et al. 119831 Pethick & Ravenhall L1995J adjacent to the 
stellar core. The mantle is predicted only by some mod- 
els of nucleon-nucleon interaction (and is not predicted by 
other models). Thus, the existence of the mantle is hypo- 
thetical. If exists, it cannot noticeably affect the equation 
of state, and the hydrostatic NS structure (particularly, 
NS masses and radii). We expect that the strongest man- 
ifestation of the mantle consists in opening direct Urea 
process. It can operate in the two last phases of nonspher- 
ical nuclei (inverted cylinders and inverted spheres), where 
the continuum proton spectrum is formed (e.g., Oyamatsu 
1993J. The emissivity of the new process in a nonsuper- 
fluid mantle appears to be 2-3 orders of magnitude higher 
than the neutrino emissivity in the nonsuperfluid outer 
NS core. 

We have performed illustrative calculations of NS cool- 
ing which show that the new process can noticeably affect 
the cooling of low-mass NSs. Its effect is most pronounced 
in NSs with strongly superfluid cores (to reduce the neu- 
trino emission from the cores) and nonsuperfluid mantles 
(to fully open direct Urea process there). Thus, direct Urea 
process in the mantle represents a new regulator of the 
cooling of low-mass NSs. 

Our calculation of the emissivity of the new process 
is simplified (based on the Thomas-Fermi model with a 
simplified form of scalar nucleon interaction, and an ap- 
proximate choice of nucleon Fermi momenta) . The calcula- 
tion can be improved but we expect that the main results 
will be qualitatively the same. One cannot exclude (Jones 
2001) that direct Urea process operates also in the crust 
of spherical atomic nuclei, or in some selected layers of the 
crust, but its calculation is difficult (requires exact wave 
functions of nucleons). If operates, it could be a stronger 
regulator of NS cooling than direct Urea process in the 
mantle. 

It is important that delicate properties of subnuclear 
matter can potentially be tested by observations of cooling 
NSs. As clear from our discussion, NSs hottest for their 
ages are the most suitable targets of such tests. 
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